Influence of boundaries on pattern selection in through-flow * 



G^ 
(N 



D. Roth, P. Biichel, M. Liicke, and H. W. Miiller 

Institut fiir Theoretische Physik, Universitdt des Saarlandes, D-66O4I Saarbriicken, Germany 

M. Kamps 

Hdchstleistungsrechenzentrum, Forschungszentrum, D-52425 Julich, Germany 

R. Schmitz 

Institut fiir Festkorperforschung, Forschungszentrum, D-52425 Jiilich, Germany 

The problem of pattern selection in absolutely unstable open flow systems is investigated by 
considering the example of Rayleigh-Benard convection. The spatiotemporal structure of convec- 
tion rolls propagating downstream in an externally imposed flow is determined for six different 
inlet/outlet boundary conditions. Results are obtained by numerical simulations of the Navier- 
Stokes equations and by comparison with the corresponding Ginzburg-Landau amplitude equation. 
A unique selection process is observed being a function of the control parameters and the boundary 
conditions but independent of the history and the system length. The problem can be formulated in 
terms of a nonlinear eigen/boundary value problem where the frequency of the propagating pattern 
is the eigenvalue. 
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I. INTRODUCTION 



The spontaneous formation of patterns in dissipative 
systems with a continuous energy supply is a commonly 
observed phenomenon |^]. Classical examples are the 
generation of rolls in Rayleigh-Benard convection (RBC) 
and Taylor vortex flow (TVF) in the annulus be- 
tween two concentric cylinders [Q,^ . If the forcing is suf- 
ficiently strong patterns bifurcate out of a homogeneous 
base state and often show a periodic structure in space 
and/or in time. For a given set of parameters and bound- 
ary conditions (BCs) one frequently observes a discrete 
or continuous spectrum of allowed wave numbers and fre- 
quencies. The initial conditions finally determine which 
ones are selected. Usually secondary instabilities limit 
the range of stable existence in parameter space. Recent 
numerical simulations of RBC [|-§ and TVF suggest 
that this degeneracy of stable stationary patterns is lifted 
as soon as a flow through the system is imposed. These 
"open flow systems" exhibit a pattern selection mecha- 
nism which uniquely fixes the dissipative structure. 

Spatially localized perturbations in open flow systems 
are advected by the flow. Depending on whether the 
trailing front of a growing perturbation of the basic state 
propagates within a linear analysis in the downstream 
or in the upstream direction the system is convectively 
or absolutely unstable, respectively The property of 
absolute or convective instability can be controlled exter- 
nally by tuning the strength of the imposed flow [p|-p2| . 
However, this distinction was made only in the more re- 



cent experiments [ pd]jl^ ]. Pattern formation in convec- 
tively unstable systems requires a permanent source of 
disturbances in order to sustain steady convection [ p^ . 
The spatiotemporal behavior of such "noise sustained 
structures" sensitively depends on the details of the per- 
turbation source [ pdUl^ . 

The present paper, however, deals with the absolutely 
unstable situation where the nonlinear structures are self 
sustained and stable. They are characterized in the 
longtime limit by time periodic fields with characteristic 
streamwise field profiles and wave number profiles. The 
frequency associated with the downstream motion of the 
pattern turns out to be constant in space and time. For 
a given set of control parameters and BCs, the pattern in 
the absolutely unstable regime is independent of initial 
conditions and history. 

In this paper we deal with RBC subject to an imposed 
lateral shear flow, i.e., with convection rolls propagat- 
ing in downstream direction. A related investigation of 
TVF has been published elsewhere 0| . Different kinds of 
inlet and outlet BCs are considered and their influence 
upon the convection structures is systematically inves- 
tigated. To this end we perform computer simulations 
of the full hydrodynamical field equations that we re- 
fer to as Navier-Stokes equations (NSE) for short. The 
results are compared with analytical and numerical so- 
lutions of the corresponding Ginzburg-Landau equation 
(CLE). The amplitude equation most clearly reveals that 
the selection process can be understood as a nonlinear 
eigenvalue problem ||^: The frequency plays the role of 
the eigenvalue, while the set of hydrodynamic fields is the 
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associated eigenfunction. Then the selection seems to re- 
sult from requiring the spatial variation of pattern enve- 
lope and phase to be as small as possible under the im- 
posed inlet/outlet BC which is analogous to ground state 
properties of the linear stationary Schrodinger equation. 

The paper is organized as follows: Section II speci- 
fies the system, the BCs for the simulations, and the 
methods of investigation. Section III recapitulates pat- 
tern selection properties in the presence of through-flow 
within the GLE. This provides a useful frame work for 
the interpretation of our results. The latter ones are pre- 
sented in Section IV which forms the main part of our 
paper. A comprehensive analysis of the convection pat- 
terns is given and a discussion of the influence of BCs 
on spatiotemporal properties of the selected structures 
as a function of through-flow rate. The conclusion gives 
a short summary of our results. 



II. THE SYSTEM 

We investigate properties of convection rolls in a vis- 
cous, incompressible fluid layer, that is heated from below 
in a homogeneous gravitational field directed downwards. 
In addition we impose a lateral through-flow. 



measures the heating and the Reynolds number 



Ud 

V 



(2.5) 



A. Equations 



measures the through-flow rate. Here a denotes the ther- 
mal expansion coefficient of the fluid, g is the gravita- 
tional acceleration, and U the vertical average of the im- 
posed lateral through-flow. 



B. Conductive state 

For small Rayleigh numbers one observes a homoge- 
neous basic state with a linear conductive temperature 
profile Tcond and a parabolic Poiseuille profile lJ{z) of 
the velocity field 

■o-cond ^ U {z)^x ; U {£) = ^aRez(\ ~ z) ; (2.6) 
Tco„d = --z, (2.7) 

which is globally stable for small thermal stress. Note 
that the effective through-flow strength is given by the 
Peclet number Pe = aRe. We have investigated in this 
work the case ct = 1. Thus, the comparison of through- 
flow induced effects in fluids of other cr has to be based 
on common Peclet numbers. 



The system is described by the balance equations for 
mass, heat, and momentum in Oberbeck-Boussinesq ap- 
proximation [0|4|-[l|] 

= -V-u, (2.1) 



dtT 



-V Q ; Q = uT-VT, 



dtVL= —V (u : u + p — (T V 



(2.2) 



(2.3) 



that we shall refer to as Navier Stokes equations (NSE). 
Lengths are scaled by the thickness of the layer d, times 
by the vertical thermal diffusion time (P / k, and the ve- 
locity field u — (u, V, w) by n/d. Here T denotes the 
deviation of the temperature from the mean tempera- 
ture To in the fluid measured in units of the temperature 
difference AT between the plates and p is the reduced 
pressure. The ratio of the momentum diffusivity ly and 
thermal diffusivity k, i. e., the Prandtl number a = v/k 
is the material parameter characterizing the fluid. The 
through-flow is in positive x-direction, with the inlet at 
X = and the outlet a,t x — L. In the vertical direction 
the lower plate is located at z — 0, whereas the upper 
one is at z = 1. 

The system is characterized by two control parameters: 
The Rayleigh number 



Ra = 



agd^AT 



(2.4) 



C. Transverse convection rolls 

At a critical Rayleigh number RadRe) |lj,^,|| a lat- 
erally periodic solution describing downstream traveling 
transverse convection rolls with axes perpendicular to the 
through-flow direction bifurcates out of the basic state. 
For convenience we introduce the relative deviation 



Ra 



Rac (Re) 



1 



(2.8) 



from the threshold as a new parameter measuring the ex- 
ternally imposed thermal stress. We also use the reduced 
distance 



Ra , 

e — K — 1 



(2.9) 



from the threshold Rar^ for zero through-flow such that 



M = 



1 + Cc (i?e) 



(2.10) 



with ec = Rac{Re)/Ral - 1. 

In this work we will focus on patterns in which the 
roll axes are perpendicular to the through-flow direc- 
tion. They occur as the first instability in narrow 
convection channels which suppress longitudinal rolls 
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9 10 2(1,G1| whereas in laterally unbounded systems 



6u{x, z\ t) = u{x, z\ t) — U{z) 



(2.14) 



the bifurcation threshold for longitudinal rolls is lowest 
pTf . In our numerical simulations we investigate the x-z 
plane perpendicular to the roll axes ignoring any spa- 
tial variation along the roll axes. Thus we solve the 
hydrodynamical field equations (2j]-2T3|) for the lateral 
velocity component u{x,z]t), the vertical velocity com- 
ponent w{x, z; t), the pressure p(x, z\ t), and the temper- 
ature field T(x, z; t). 



D. Boundary conditions 

The horizontal BCs are always no slip and perfectly 
heat conducting 

u = w = atz = 0, 1 (2.11) 
T = 1/2 atz = 0, T = -1/2 atz = 1. (2.12) 

At the lateral boundaries we always enforce the Poiseuille 
flow profile 



u{x, z;t) = U(z) at x = 0,L. 



(2.13) 



For the vertical velocity we have applied laterally no 
slip (N) conditions (cf. Table |) or free slip (F) condi- 
tions. For the temperature field we have investigated 
the three lateral conditions that are listed in Table || and 
that are abbreviated by the symbols Tcond, Tq, and Qo, 
respectively. We have called the latter thermally insulat- 
ing since the diffusive part of the heat current Q (2.2) 
vanishes there. But it should be noted that the lateral 
Poiseuille flow transports heat advectively through the 
lateral boundaries. To specify the six possible combi- 
nations of {N,F} conditions with {Tcond,To,Qo} con- 
ditions, that we have investigated, we use a two-letter 
coding. Thus, e. g., N-Qo implies w — and dxT — 0. 
All numerical solutions have been obtained for a channel 
length of L = 50. 



E. Methods of investigation 



To solve the time dependent 2D NSE ( |2T[]2^ ) we used 
an explicit finite difference numerical code that allows 
a flexible and convenient incorporation of various BCs 
and that has been used sucessfuUy in related problems 
. The final-state solutions obtained in this way rep- 
resent transverse convection rolls that propagate down- 
stream under a stationary x-dependent intensity enve- 
lope. These patterns oscillate in time with a global, i. e., 
spatially independent frequency u while wave number k, 
phase velocity Vph = uj/k , and convection intensity vary 
in a;-direction. The frequency and the spatial variation 
of the final pattern depend on Ra, Re, and the BCs but 
not on parameter history. 

To characterize the convection structures we use the 
convective deviation 



of the lateral velocity field from the Poiseuille flow as the 
order parameter field. Compared with the vertical veloc- 
ity w or the convective temperature profile T — Tcond this 
field has some technical advantages (cf. Sec. [V). After 
having determined the global oscillation frequency uj of 
the pattern we perform a temporal Fourier analysis 



Su{x, z; t) = Sun{x, z)e 



(2.15) 



The zeroth component 5uq{x,z) describes a stationary 
modification of the basic state's Poiseuille through-flow 
profile that reflects the presence of a stationary secondary 
roll vortex system (cf. Sec. iV). The higher harmonics 
Un with n>2 are in general small compared to the fun- 
damental mode 



ui{x,z) = |Mi(a;,z)|e*'^i(^'^). 



(2.16) 



We use |mi I to characterize the intensity profile of the pat- 
tern and we characterize its wave nu mbe r by fc = dx'fi- 
Within the CLE approximation (Sec. Ill) the fields oscil- 
late harmonically in time and only |wi| is nonzero. Thus 
the CLE cannot reproduce the stationary secondary vor- 
tex flow discussed in Sec. IV. 



III. PATTERN SELECTION WITHIN THE CLE 

The selection process for convection patterns in the 
presence of a lateral flow is a nonlinear eigenvalue / 
boundary value problem. This is easiest to understand 
within the CLE approximation However, one should 
keep in mind that CLE and NSE results differ systemat- 
ically. 



A. GLE for propagating patterns 

Within the amplitude equation approximation the con- 
vection fields of the propagating pattern of transverse 
rolls have the form of harmonic waves, e. g., 

5U (X, Z-t)=A {X, t) e*('=c:'=-^ct)^ _^ ^ ^ (3 ;l) 

The complex amplitude A{x, t) is common to all fields 
and obeys the GLE 



To [A + VgA') =n{l + icq) a + ^0 (1 + ici) A" 



~j{l + iC2)\A\^A (3.2) 
subject to the homogeneous BC 

A{x,t) =0 at X = 0,L. (3.3) 
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Dot and primes denote temporal and spatial derivatives. 
This BC effectively imposes on the fields in GLE approx- 
imation the N-Tcond condition: T — Tcnnd, Su — 0, and 
w = 0. The other BCs listed in Table || are not realiz- 
able within the GLE approximation since it enforces one 
common amplitude for all fields. However, we found that 
the NSE results for envelope profiles and growth lengths 
obtained with the F-Tcond BC compare better with the 
GLE results. The GLE ansatz for the fields implies via 
the critical functions a spatially constant phase difference 
of 7r/2 between 5u and w. This is not the case for the 
NSE fields subject to the N-Tcond BC, which show close 
to the boundaries a deviation from the above phase rela- 
tion, while the F-Tcond BC entails phase relations closer 
to those of the GLE. 

All coefficients of the GLE, the critical numbers ec, 
kc^ LOc, and the eigenfunction u{z) have been calculated 
1^ as functions of Re. As a consequence of the sys- 
tem's invariance under the combined symmetry opera- 
tion {x —X, Re —Re} kc, tq, 'Coi 7 ^'^s even in Re 



and the imaginary parts 



while ujc, the group velocity 

co) ci, C2 are odd in Re P,|22||. For cr = 1 an expansion in 
Re yields « (i?e/35.93)Vfcc ~ 3.116[1 -I- (i?e/99.97)2], 
3.6475i?e, ro « 0.07693fl - (i?e/39.13)2], « 



0.148[1 - (i?e/35.24)2 



1.229i?e, 7 w 0.7027[1 



(i?e/332.6)2], Co « i?e/140.1, ci « i?e/39.98, C2 ~ 
Re/ 387 @. " The next order terms are in each case two 
orders in Re higher. 



B. Eigenvalue problem 

Let us consider small through-flow rates Re and ther- 
mal stresses ^ slightly above the boundary 



(3.4) 



between the absolutely and convectiyely unstable param- 
eter regime. We solved the GLE (|3.2| ) numerically with a 
Cranck-Nicholson algorithm with a resolution of 20 grid 
poin ts per unit length d . Th e fin al state solution of 
(|3.2D with the homogeneous BC (3.3) is of the form 



A{x,t) = a(a;)e-*"* = i?(x)e*['^(^)-"*l . 



It oscillates harmonically with the global frequency 

n = LU — LUc (3-6) 

under a stationary envelope R{x) with a stationary wave 
number profile 



q{x) — k(x) — kc = f'{x). 



(3.7) 



Inserting the solution ansatz (^^) into the GLE one ob- 
tains the nonlinear eigenvalue problem 



iro(— n -I- Vgq)R + ToVgR' = ^{1 + ic^jR 
+eo(l + lCi){R" -q^R + iq'R + 2iqR') 
-7(l+iC2)i?^ 



(3.8) 



The eigenvalue fl and the associated eigenfunction 
a{x) = i?(x)e*'^'^^-' describing dynamics and structure of 
the pattern are fixed by the BC (p.3|). We cannot make 
statements about th e eig envalue spectrum. But our nu- 
merical solutions of (|3.2| ) suggest that initial conditions 
always evolve into a pattern with a smooth amplitude 
a{x) with smal l fr equency fl. In tha t sense the eigen- 
value problem ( |3.8| ) with the BC (3.3) resembles ground 
state problems of the Schrodinger equation - increasing 
the spatial variation of the wave function (amplitude a) 
implies higher kinetic energy (frequency Q). 



IV. PATTERN SELECTION WITHIN THE NSE 
FOR DIFFERENT BOUNDARY CONDITIONS 

In this section we investigate how structure and dy- 
namics of patterns selected in different through-flows are 
influenced by the lateral BC. We present numerical sim- 
ulations of the full fleld equations in comparison with nu- 
merical solutions of the GLE. Results are often presented 
as functions of the scaled group velocity 



(4.1) 



which scales as Re/^/Jl for small Re (cf. Sec. III). With 



this scaling the boundary between absolute and convec- 
tive instability lies at = 2 and thus allows direct com- 
parisons for different fi-Re combinations. We have varied 
the Reynolds number at a fixed e = 0.215. Thus we scan 
the Re-fj, control parameter plane according to (2.1C) 
along the path fj, = 0.215/ [1 -f ec(-Re)] that is slightly 
curved and that hits the borderline fJ^^onv ( |3-4D of the ab- 
solutely unstable regime at i?e ~ 3.8. Along this path Vg 
varies for a ^ 1 as Vg 0.53i?e + 0[Re^). 



A. Selected frequency eigenvalues 



In Fig. |l| we show NSE frequency eigenvalues (symbols) 



(3.5) as functions of the scaled group velocity Vg ( [4.1[ ) for the 
six BCs of Table I together with the GLE eigenvalue (full 
line) for N-Tcond conditions. Approaching the border to 
the convective instability at Vg — 2 the growth length of 
the patterns from the inlet increases and the dependence 
of Lj on the BC gets weak and all uj values come close 
together, whereas for small Vg much stronger differences 
can be observed. The selected frequencies are predomi- 
nantly influenced by the inlet BC. A change of the outlet 
BC leads to w-variations which are 10"'^ to 10~^ times 
smaller than those induced by changing the inlet BC. 
This insensitivity of w to the outlet BC reflects the differ- 
ent amplifying properties of perturbations in upstream or 
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downstream direction. However, the intensity envelopes 
(Sec. IVB) and the wave number profiles (Sec. ^v|d) are 
more strongly influenced by the outlet BC, in particular 
in its vicinity. 

From all BCs N-Tcond leads to the largest and F-Qo 
to the smallest frequencies. In general a change in the 
inlet BC from no slip to free slip causes cu to decrease. 
For Tcond and Tq conditions this change is small, whereas 
it becomes significant for the thermally insulating cases 
Qq. This can be explained by the different phase pin- 
ning properties of the different BCs. In the N-T^ond case 
that enforces the basic state at in- and outlet the tem- 
poral fourier amplitudes of all convective fields drop to 
zero, leaving the phase indetermined and free to move. 
On the other hand, enforcing a von Neumann BC for 
one of the fields yields a nonvanishing field amplitude, 
a phase pinning effect, and therefore a reduction of the 
phase velocity. In the case of F-Qq the phase pinning 
is strongest, because there dxW as well as d^T vanish, 
leaving a free-phase condition only for u. 

Both To conditions induce a large stationary "bound- 
ary vortex" flow in the vicinity of the boundaries, that 
can be seen directly via the zeroth Fourier modes of the 
fields, e. g., Suq and indirectly via the lateral variation of, 
e.g., |wi| (cf. Sec. IVB and Fig. ||). It is quite similar to 
the stationary Ekman vortex structures that appear in 
the Taylor Couette system with through-flow as a result 
of a vanishing azimuthal velocity at inlet and outlet 0. 
Besides that the lateral variation of induces for all 



BCs weaker stationary "flank vortices" (cf. Sec. |Vp3 and 
Fig. 1^) localized in the region of largest lateral growth 
of The stationary vortex structures exert an ad- 

ditional frictional force on the phase of the propagating 
vortices that reduces their phase velocity, when the fields 
Suq and ui spatially overlap. The frequency reduction 
roughly depends on the (Suo-amplitudes of the station- 
ary secondary vortex structures relative to the size of the 
propagating vortex amplitudes, sa y An analogous 

argument was given by Cross et al. |23|| for the wave num- 
ber selection at Re = 0. In this way one can understand 
(i) That the tj-reduction for Tcond and Qo boundaries 
saturates with decreasing Vg-. For these BCs the mag- 
nitude of the i5Mo-amplitudes of the "flank vortices" ap- 
proaches near the boundaries the size of from below 
when the |Mi|-profile of the propagating structure comes 
closer and closer towards the inlet for ^ 0. (ii) That 
u) decreases for small Vg dramatically for To-boundaries: 
They enforce stationary "boundary vortex" amplitudes 
Suo > \ui\ such that (5uo/|ui| further increases when 
— + 0. As an aside we mention that the decrease of 
uj/uJc—1 at small through- flow in the Tq case differs from 
the otherwise very similar results for propagating vortex 
patterns in the rotating Couette system with Ekman vor- 
tex enforcing inlet/outlet conditions. There an increase 
of k/kc — 1, i- e., of w/wc — 1 has been observed for small 
Re 0, Fig. 6]. 

The CLE frequency eigenvalues (full curve in Fig. |^) 
deviate in a systematic manner from the NSE ones (open 



circles for the common N -Tcond BC). This has also been 
found for downstream propagating vortices in the Taylor 
Couette system j?] . The reason for this systematic devi- 
ation seems to be here as for propagating Taylor vortices 
the different dispersion relations of NSE and CLE. See 
Ref. for a detailed discussion. 



B. Lateral pattern profiles 

In Fig. ^ we present the lateral variations of and 
5uq at 2 = 0.225 for two different Re values. The bulk 
region 7 < a; < 45 is cut out, because there the Fourier 
modes are almost constant. Since 5u is zero at in- and 
outlet for all our BCs this field allows best to elucidate 
the changes induced by different BCs. The lateral ve- 
locity 5u has a further advantage for the analysis of the 
BC-dependence of convection patterns since 5u is more 
sensitive to lateral BCs than, say, w: The lateral defor- 
mations in the flanks of the profiles of 5u are more pro- 
nounced than for w. In our earlier work where w 
profiles were analyzed we detected only minor differences 
between NSE and CLE profiles where these deformations 
are completely absent. 

The deformations in the lateral field profiles of, say 
|ui I, are caused by the nonlinear couplings to a stationary 
vortex structure located in the flank regions where the 
lateral variation of the intensity of the propagating 
pattern is large. Typically the secondary flow amplitude 
5uq of the "flank vortices" reaches a strength of about 
2% of the bulk The "flank vortex" flow increases 

when changing the BCs according to — > F and/or 
Tcond ^ To ^ Qo- The "flank vortices" associated with 
the inlet facing flanks of |ui | become weaker with increas- 
ing Vg since then the |Mi|-profile becomes smoother as its 
growth length increases (see Fig. ||d for Re = 8/3). 

Besides this |mi [-induced stationary "flank vortex" 
flow, there is for the Tq BC a boundary induced sta- 
tionary "boundary vortex" structure of large amplitude 
Suo ~ which is very similar to the stationary Ek- 
man vortex structure in the analogous Taylor Couette 
system [Q. This "boundary vortex" structure, induced 
by the Tq BC, is always located near the lateral bound- 
aries irrespective where the flanks are located and its 
intensity falls off towards the bulk - see Suq in Fig. ||d 
for Re = 8/3. 

For a more quantitative comparison of the selected am- 
plitude's profiles we have determined the characteristic 
growth lengths from inlet and outlet where |wi| reaches 
half of its bulk value. To that end we have performed a 
least square fit to with the function (1 — e''^)/(6-f e°^) 
with fit parameters a and b in order to obtain also in the 
presence of amplitude deformations comparable results. 
This fit function yields good results except directly in the 
inlet region. In Fig. ^ we show the scaled growth lengths 



(4.2) 
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of as a function of Vg. Note that inlet (outlet) 
growth lengths are shown for positive (negative) Vg. L 
has to be identical at in- and outlet for Vg — 0. In ad- 
dition L seems to have there a continuous first deriva- 
tive for Tcond and Qq conditions but not for Tq. For 
the latter BC L strongly increases as Vg approaches zero 
since the stationary "boundary vortex" induced by the 
To condition strengthens with decreasing \Vg\ relative to 
I Ml I and therefore pushes the oscillating pattern into the 
bulk. Changing the BCs along the sequence F-Tq — > A^- 
To ^ iV-Tco„<i ^ F-T^ond N-Qo ^ F-Qo decreases 
L successively. For both Tcond conditions GLE and NSE 
agree quite well with deviations being smaller for F-Tcond 
BC than for N-Tcond BC (see Sec. p|A). 



which seems to be at first sight at odds with the GLE 
formula Rf — ^ — ^^g^ is caused by two facts: (i) Due to 
the pattern selection process the selected wave number 
in the bulk region is uniquely determined by the control 
parameters /i, Re. For a give n /i one obtains for Re ~ 
the critical wave number |23|| . Then, with increasing Re 
the wave number first slightly increases and later on it 
decreases below the critical one to the value 



^0% 



(4.5) 



C. Phase relations in the bulk 



Cl - C2 

at the border between absolute and convective instabil- 
ity 0- (ii) The contribution from ^gg^ to i?^ is very 
small since \qh\ ^ 0.05. Therefore the variation of i?^ 
and thus of is mainly caused by monotonously de- 
creasing [I — 0.215/[1 -f (i?e/35.93)^] while increasing Re 
from Re — Q to Re 3.8 where Vg = 2. 



Let us consider first the dispersion relation between 
the frequency w, i.e., the temporal phase gradient and 
the spatial gradient fc in the bulk part of the system. 
The GLE (1^) predicts the relation 



(4.3) 



Tofi = (C2 - Co)^* + ToWg(7b + (Ci - 02)^0*; 



between selected frequency £7 = 0; — cjc, wave number 
Qb = kb — kc, and modulus i?^ = (/i — Co'?6)/7 ^'^^ bulk 
patterns with = R'^ = = 0. Fig. |a shows the GLE 
result ( L^ ) (full line) together with the NSE results (sym- 
bols) obtained for the six BCs for all Reynolds numbers 
sufficiently large to establish a bulk region with spatially 
constant wave number. Since the bulk GLE wave num- 
bers are very close to the critical one the quadratic con- 
tribution to fl and to Ri, is for all through-flow rates 
< < 2 very small so that O is practically linear in 
k (Fig. |a). That also holds for the NSE case. Here, 
however, the spread in the bulk wave numbers around kc 
is considerable larger depending on the lateral BC and 
Re. Within the resolution of Fig. 4a the bulk dispersion 
relation for propagating roll vortices resulting from GLE 
or NSE reads 



LOc + Vg{k ~ kc) 



(4.4) 



irrespective of the BC. The slight deviation between the 
thin dotted curve in Fig. ^a that represents uj/uJc— 1 — 
(fc— kc)vg/u!c according to (4.4) and the GLE result (full 
line) comes from the contrib utio n (c2 — Co) ^ / {lh cTq) to 
the GLE dispersion relation ([t.3|). 

In Fig. 4b we show the relation between bulk intensity 
of the propagating pattern and bulk wave number. 
Again the symbols (full hue) denote the NSE (GLE) re- 
sults for several Re a nd B Cs. To determine jwip for the 
GLE we have used ( |3.l|) with the properly normalized 
linear eigenfunction u{z) |24j evaluated at Re — 0. The 
relations between bulk pattern amplitude and bulk wave 
number resulting from NSE and GLE are obviously quite 
different. The peculiar shape of the full curve in Fig. ||b 



D. Wave number profiles 

Having investigated the frequency eigenvalue and the 
modulus profile |wi| of the associated eigenfunction we 
now address its phase structure. The frequency eigen- 
value and the lateral profile of the phase gradient k{x) = 
dx(pi{x) seem to be closely tied together. And the latter 
is largely slaved by the envelope profile which itself 
is most strongly influenced by the BCs. 

Fig. 1^ shows selected modulus profiles of |ui| together 
with corresponding wave number profiles k{x) for a small 
through-flow with Re — 1/6. The inlet boundary condi- 
tions were N-Tcond or F-Qq and the outlet conditions 
were always N-Tcond- The frequencies for these two inlet 
BCs differ significantly (c.f. Fig. 1) and also the wave 
number profiles k(x) are quite different. Large deforma- 
tions in the pattern amplitude cause strong deforma- 
tions of the phase gradient k(x). So both, the temporal 
derivative to and the spatial gradient k{x) of the phase of 
the propagating pattern are influenced by the inlet BC. 
Note also that the (boundary induced) deformations of 
the pattern profile correlate with corresponding vari- 
ations of k{x). The GLE, on the other hand, yields a 
wave number profile (open circles in Fig. 5b) that is very 
close to kc and that is practically constant (c.f. Fig. 4b] 
for more details) in space while the NSE wave numbers 
show significantly larger variations. 

In the inset of Fig. 5 we show k{x) in the bulk region 
10 < a; < 45 where is almost constant. For the phase 
pinning F-Qq BC the wave number does not saturate 
into a bulk part for the small through-flow of Re = 1/6. 
Instead k(x) shows an exponential growth (full line in the 
inset of Fig. 5) in the region of bulk pattern amplitude 
from small k in the upstream portion to large k in the 
downstream region: For sufficiently small through-fiow 
the phase pinning forces at the inlet are strong enough 
to cause an elongation of the pattern, i.e., a decrease in 
k near the inlet that heals out exponentially in down- 
stream direction. Increasing Re or decreasing the inlet 
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phase pinning with other BCs reduces the pattern elon- 
gation, i.e., the wave number gradients. On the other 
hand, decreasing Re increases the gradient of k{x) and 
eventually even suppresses the phase propagation alto- 
gether. This latter transition at very small Re depends 
on the system size — c.f. Rcf. 1^ for a related detailed 
investigation. 

However, when the parameters are such that a true 
bulk part of the pattern is established in which ampli- 
tude and wave number are constant then the selected 
frequency and lateral pattern profile seem to be indepen- 
dent of the system size. This always occurs in sufficiently 
large systems when the through-flow is sufficiently strong. 



V. CONCLUSION 

We have investigated convection roll structures in sys- 
tems of finite length that propagate downstream in an ex- 
ternally applied lateral flow with roll axes perpendicular 
to the flow direction subject to six different inlet /outlet 
BCs. Within the absolutely unstable parameter regime 
a unique pattern selection is observed. The selected con- 
vection patterns are independent of parameter history, 
initial condition, and system length provided the latter 
is large enough to allow for a saturated bulk region with 
homogeneous pattern amplitude and wave number. But 
they depend on the lateral BC. The outlet BC influ- 
ences the convection structures only locally in the vicin- 
ity of the outlet while the globally constant frequency 
and the lateral profiles of pattern amplitude and wave 
number and their bulk values are basically determined 
by the inlet BC. The selected frequency uj of the pat- 
tern oscillation is an eigenvalue of a nonlinear eigenvalue 
problem with the set of time periodic convection fields, 
/(x, z, t) = f{x, z,t + 27r/w), being the associated eigen- 
function. The analysis of the appropriate CLE approxi- 
mation suggests that the eigenfunction associated to the 
selected eigenfrequency varies as smoothly and as little 
as possible under the imposed BC. This property is sim- 
ilar to the ground state behavior of the linear stationary 
quantum mechanical Schrodinger equation. 

Different lateral BCs entail different frequencies and 
eigenfunctions, i.e., lateral pattern profiles. For example, 
inlet conditions that exert a phase pinning force on one or 
more of the fields of traveling convection rolls reduce the 
frequency. The spatial growth behavior and the pattern 
profile of the convective structures in the region between 
inlet and bulk saturated pattern is strongly influenced 
by the inlet BC and can be quite intricate. While the 
frequencies and bulk wave numbers obtained for various 
BCs and through-flow rates well inside the absolutely un- 
stable regime differ substantially they all fall onto a lin- 
ear dispersion curve. Finally, when approaching with in- 
creasing Re the border between absolute and convective 
instability the eigenfrequencies and also the spatial pat- 
tern profiles belonging to different BCs approach each 



other — with increasing growth length of the pattern 
from the inlet the influence of the latter on spatiotempo- 
ral properties of the convection pattern becomes weaker. 
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TABLE I. Lateral boundary conditions at a; = and/or 
X = L. 



field 


type 


symbol 


w{x, z;t) = 


no slip 


N 


dxw{x, z;t) = 


free slip 


F 


T{x,Z;t) = Tcond{z) 


conductive profile 


Tcond 


T{x,z;t) = 


mean temperature 


To 


d^T{x,z;t) = 


insulating 


Qo 



FIG. 1. Selected oscillation frequency of traveling convec- 
tion rolls obtained from the NSE for different lateral BGs 
(symbols) and from the GLE fo r N -Tcond BC (full line) ver- 
sus the scaled group velocity Vg ( [4.l[ ). Dashed lines are guides 
to the eye. Parameters are e = 0.215 and a — 1. 

FIG. 2. Lateral profiles of convection structures for differ- 
ent BGs. Temporal Fouriermodes |wi| and Suq reduced by the 
bulk value of |ui| at x = 30 are shown ai z — 0.225 close to 
the inlet {x — 0) and outlet {x = 50) boundaries. For clarity 
Suq is suppressed in the outlet region 47 < a; < 50. Lines 
refer to NSE results and open circles to the GLE subject to 
N-Tcond conditions. Parameters are e = 0.215 and a = 1. 

FIG. 3. Scaled growth lengths L of propagating convection 
patterns versus scaled group velocity Vg. Growth distances 
where |wi| reaches half its bulk value from inlet (outlet) are 
shown for positive (negative) Vg. Symbols refer to NSE results 
and line to GLE. Parameters are e — 0.215 and (7 = 1. 

FIG. 4. (a) Frequency eigenvalues lo/uc — 1 and (b) 
bulk flow intensities versus selected bulk wave num- 

bers — fee of propagating roll patterns are shown for sev- 
eral Reynolds numbers covering the whole absolutely unstable 
regime < Vg < 2. Symbols (lines) refer to results obtained 
from the NSE with different BGs (GLE with N-Tcoud BG). 
The thin dotted curve in (a) shows {k — kc)vg/uc as discussed 
in Sec. IV G. Parameters are e = 0.215 and a = 1. 

FIG. 5. (a) Amplitude profiles and (b) wave number pro- 
files close to inlet and outlet. The temporal Fourier mode |mi| 
of the traveling convection pattern was reduced by the bulk 
value at a; = 30. The inset shows the lateral variation of the 
wave number in the center region of the system. The outlet 
BG was N-Tcond, the inlet BGs are indicated in the inset. 
Parameters are e = 0.215 and a = 1. 
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Influence of boundaries on pattern selection 
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Influence of boundaries on pattern selection 
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Influence of boundaries on pattern selection 



D. Roth et al., Physica D 



